Effect of H ion formation on Positronium-Hydrogen 

elastic scattering 

P K Biswas 

Departamento de Ffsica, Instituto Tecnologico da Aeronautica, CTA 
Sao Jose dos Campos 12228-900, SP, Brasil 
February 2, 2008 

Abstract 

Effect of charge-transfer recombination for positronium (Ps) scattering is studied 
on Ps-H system using coupled-channel formalism considering a new rearrangement 
channel Ps+H— > e + +H~ and exchange. The correlation and continuum effects in- 
troduced by this charge-transfer channel result to substantial convergence to the low 
energy scattering parameters. Effects on scattering length, PsH binding energy, and 
low-energy (0-10 eV) cross sections are evaluated and H~ formation cross sections 
are reported from above the threshold (6.438 eV) to 100 eV. 
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Introduction : Charge-transfer recombinations have been found to be of fundamental 
importance in simulating the continuum effect in the scattering dynamics of positrons 
[H, |2j. Their roles are yet to be investigated for the scattering of positronium (Ps) atom - 
an exotic bound state of electron and positron with singlet and triplet spins. This seems 
to be very important since, despite a long history of theoretical studies |3|, |J || || 0, ||, 
HH m, 0, m, 0, |3, P H, H, ||, ||, ||, we are yet to produce 



converged scattering results in the ab initio coupled-channel (CC) framework, even for 
the simplest Ps-H system [TOI . 



So far, in dealing with Ps scattering, emphasis has been given to the exchange rear- 
rangement of electrons, which is also of fundamental importance, as the static-potential 
for a Ps-target system vanishes due to internal charge-and-mass symmetry of Ps. However, 
exchange constitutes only a part of the continuum and Ps exhibits a strong polarizability 
of 36ag, which is expected to make it vulnerable to continuum effects and charge transfer 
processes. Continuum effect can be simulated through pseudostates of Ps and target, 
but for Ps scattering this is expected lead to an untractable situation as the number of 
scattering channels grows like iV 2 , where A^-number of intermediate states for each of Ps 
and target are considered. Considering the difficulty, so far, the pseudostate technique 
has been applied (to Ps-H and Ps-He scattering problem) employing pseudostates for only 
Ps. This is not expected to simulate the appropriate continuum effect as the target is 
kept idle to keep the calculational schemes tractable. However, from the experience in 
positron-atom scattering g], we understand that allowed charge transfer processes can 
also enhance the continuum effect significantly and lead to the required convergence of 
the theoretical scheme. The effect of such charge-transfer channels is yet unexplored for 
Ps scattering and is the subject of this work. 

Apart from the initial works B [|, || [71 0], the recent theoretical development on 
Ps scattering studies, has contributions from the Belfast group |], [TO, O, the Calcutta 



group and the Sao Paulo group |7[g§^gg,g5^^^g§gg. While 

the Belfast and Calcutta groups use complete ab initio coupled-channel (CC) formalism, 
the Sao Paulo group indulges to model exchange correlation potentials (to enhance the 
correlation and continuum effects) in the ab initio CC formalism. In this investigation, we 
concentrate on the simplest but the richest Ps-H scattering problem, where predictions 
of the ab initio coupled-channel (CC) calculations are yet to converge with those of the 
accrrate variational calculations [[!(], |28|j . 

In the ab initio CC formalism, the largest continuum effect has been taken in a 22- 
state calculation which uses three eigen states and nineteen pseudostates of Ps keeping 
the target idle. This calculation yields PsH binding energy of 0.634 eV and S-wave singlet 
resonance energy of 4.55 eV compared to the most accurate variational estimates of 1.067 



eV |30| and 4.003 eV p8[ , respectively. The lack of convergence was expected as the virtual 
effects of the target are neglected. Interestingly, predictions of this CC calculation come 
to close agreement with the predictions of another variational calculation [0] that has been 
performed without explicit electron-electron correlation in their trial ansatz. The precise 



variational calculations above |30], £8[ includes explicit e — e correlation in their trial wave 
functions and they also contain H~ flavour. These clearly emphasize the necessity of the 
introduction of more electron-electron correlation and continuum effects to the present 
state of the art CC calculations. We understand that, for Ps-H system, electron-electron 
correlation and continuum effects could be enhanced through the formation of H~-ion via 
Ps+H— > e + +H~ channel in the CC model. 

Before finishing the introduction, we note another relevant developement in this sub- 
ject. For the last few years we have been engaged in addressing the Ps scattering problems 
by using a model exchange correlation potential in the ab initio CC formalism. Apart 
from the recent theoretical development on Ps scattering by the Belfast and Calcutta 
groups employing complete ab initio formalism, studies has contributions from two more 



groups other than the Belfast group which has considered the Ps-pseudostate technique 
[P, [T(J O, mentioned above. The Calcutta group has addressed the Ps-H problem in var- 



ious ways JL2], [13], [T^| , but they have not reported the binding or resonance parameters. 
So, these aspects and hence the overall convergence of their low energy scattering results 
could not be assesed comparing with the accurate variational predictions. 

In the last few years, we are addressing Ps scattering problems by using a model reg- 
ularized exchange-correlation potential (non-local) PI] in the exact coupled Lippmann- 



Schwinger equation |7|, [TJ, [TJ, |2D|, f^, f£| fZ§ f£| [2§. In these calculations, con- 
vergence and agreement with accurate variational predictions are achieved by employing 
a minimum number of effective scattering channels and tuning the exchange correlation 
potentials by means of a parameter. The model agreed well with the measured data 
|32| , |33| , |34|| in Ps-He []. When applied to Ps-H, it reproduces the precise binding and 
resonances energies. Thereafter, it has laso been applied to other taegets (H 2 , Ne, Ar, 
Li) and agreement with measured data is obtained. However, the crucial aspect of the 
model remains the tuning of the parameter. It appears that the exchange correlation 
potentials, when tuned, compensate for the correlation and continuum effects (see details 
in the results and discussions). Thus, we expect that the direct inclusion of correlation 
and continuum effects in the model CC framework, through the H~ formation channel, 
would diminish the role of the parameter. Also, the charge transfer channel, in general, 
is expected to improve the results obtained in the ab-initio framework due to possible 
enhancement of correlation and continuum effects. 

To have a preliminary assessment of these assertions, we re-investigate the Ps-H scat- 
tering including the process Ps+H-^ e + +H~. In the present work, for a first estimation, 
we use a simple wave function of Chandrasekhar for the H~ ion . Any elaborate wave 
function containing the term would certainly enhance the electron-electron correlation 
further. We estimate the effect of this charge-transfer recombination channel in the ab 



initio CC framework by taking exchange exactly and also by using the regularized model 
exchange [^1|] in the CC equations. 

Theory: Antisymmetrizing for the constituent electrons, we expand the total wave 
function of the Ps-H system in terms of Ps (xv), H (</>/J, and H~ (if) p ) states as 



« s (n,r ai x) = EE 



+ <W? p (x)^ p (ri,r 2 ), (1) 

where r 1; r 2 denote the electron coordinates, and x is the positron coordinate of Ps; 
Si = (i"i + x)/2, ti = (i*i — x); i = 1, 2. S is the total electron spin of the system in a 
particular channel which can have values and 1, corresponding to singlet and triplet 
scattering; 5$o, is the Kronecker delta. F Vfi is the continuum orbital of Ps and G p is that 
of the rearranged positron with respect to the center of mass fixed in the target nucleus. 
The total Hamiltonian of the system may be taken as 

H = -\v 2 Ps + H° Ps + H° H + V(l) (2) 

or 

H=-~Vl + H° H -+V(2) (3) 

where H Ps , Hfj, and are the unperturbed Hamiltonian for the Ps, H, and H , 

respectively; — \V Ps and — |Vp represent the kinetic energies of the Ps and the positron, 
respectively. V(l) is the interaction potential in channel-1 (elastic) and V(2) is the same 
in the charge-transfer rearrangement channel (channel-2: e + H~). For channel-1, we have 
direct scattering (denoted by superfix d) and exchange scattering (denoted by superfix e) 
and for the transfer channel, there is no notion of exchange and we have only the direct 
scattering. The coulomb potentials for these channels are given by: 



v e d) = ! - _ - i + — (5) 

X |X — Ti| T2 |ri — r 2 | 

V(2) = (6) 
x |x — r 2 | |x — ri| 

We use exact wave function for Ps and H, and for H~, we use the wave function of 



Chandrasekhar [35]. Understanding the difficulty in dealing with a coulomb wave for 
the outgoing positron, we treat it with a plane wave in the input Born matrix element 
for channel-2, as the distortion to this plane wave will partly be done through the use 
of dynamical Lippmann-Schwinger equation (given in the next paragraph) and coupling 
over intermediate states. 

Projecting over final states of Ps (Xu 1 ), H (</y), and B.~(ip' p ), the LS equation for 
particular electronic spin state l S\ (which can have values and 1, in general but for 
e + H~ channel, it is ought to be only as H~ can exist only in singlet state), can be recast 
as IBB: 



v" n" 

r 1 f „ //^fvV(kf ; k 2 ")/^'7° M (k2", k ; 
~ ^°2^i 2 B, - k^ + *0 



f // ffiy,i/y (kf , ki // )/^/ /J // ; ^ t (ki // , kj) 



(7) 



/^^(kf,^) = ^ M (k f ,kO 

1 /•^ // jBg^(k f ,k a w )/gr° i (k a w ,k i ) 
2**$fJ dk2 fcj, _ ^ + iO 

1 \r^\r^ /* //^p'X'M"(kf> k2 // )/w"7?',i/u(k2 // > ki) 
- 2^^J dkl fcJ,-V 2 + iO (8) 

where, fc^//„« = ^r-{E — t v n — T^//}; and A; 2 ;/ = ^-{-E 1 — ty'}; and k p » are the 

on-shell momenta; m\ and m2 are the masses of e + and Ps, respectively. E represents the 
total energy of the system; e v >i, T M ", and v p " represent the binding energies of Ps, H, and 
ion, respectively. It is interesting to note that, for the charge-transfer channel, the 



summation over p" is restricted by itself as without a magnetic field, H~ can exist only in 
its ground singlet state (<S = 0) |3"7j| . We study the effect of this channel over the elastic 
scattering channel (channel-1) and thus we restrict the summation over v" and pi' to the 
ground states of Ps and H, respectively. The input potentials to the coupled equations, 
in general, are given by 

B^ M (k,,ki) = B^(kr,ki) (9) 

Where, -B d and 5 6 are the direct Born and exchange Born-Oppenheimer (BO) amplitudes. 
The exact form for B^, , , is available in the literature |21| but here we follow the sign 



convention for the Ps-wave function as used in flT8"|. For B e , we consider both the exact 



form (Appendix I) as well as the regularized form [2I| with appropriate sign convention 



[18[|. We use h = c = m 2 = 1, where m,2 is the mass of electron or positron and m\ = 2, 
the mass of Ps. 

The elastic cross section and H~ formation cross section is given by: 

a el = ^l/lsls,lslsl + (10) 
°H- = j^j^l/i°s2,l S lJ 2 (11) 

Recently, a concern has been raised by Adhikari and Mandal (AM) |38| regarding the 
validity of the calculation of the BO exchange term of Sinha et al [0 and provided a 
different results for this BO exchange term. We were also concern with AM regarding 
the feature of a minimum in the elastic cross section with BO exchange near 35 eV. 
However, present results for Bf sls lsls exchange agree exactly with those of Sinha et al 
flHf l and when this BO amplitude is used in the coupled equations, along with the required 



off-shell matrix element, it reproduces the static-exchange phase shifts of Campbell et al 
1(J, Sinha et al and Hara and Fraser Q quite precisely (discussed in the results 



and discussion). So, it implies that the BO exchange amplitudes in these calculations are 
same and the minimum in the cross section is a feature of the BO exchange term and not 
due to any error in its evaluation. AM evaluates the BO exchange term relying mostly 
on numerical calculations where the integrand contains spherical Bessel functions and 
they have reported to use 400 Gauss-quadrature points in the evaluation of the integrals. 
Whereas, using Fourier transform and Chasire integrals, we reduce the nine- dimensional 
integrals (Appendix I) to simple two-dimensional integrals which is found to converge with 
mere eight Gauss-quadrature points for each integration variable and hence any numerical 
problem is ruled out in this methodology compared to that provided by AM [38|. 



In Appendix I, we provide the analytical form of the ab initio BO exchange matrix 
element B^ sls lsls and the charge-transfer matrix element Bf g2 lsls is detailed in Appendix 
II. In a separate study, we employ the regularized model exchange for B\ sls lsls , which is 
taken from ref. [0 (it differs only in sign convention from ref. pi||). 

Numerical Procedures : In the present investigation, we consider static-exchange 
(SE) and SE plus charge-transfer reactions of Ps-H to e + -H~. The three-dimensional 
LS equations, for a particular electronic-spin state (S) are decomposed to coupled one- 
dimensional partial wave equations, which are then solved by the method of matrix in- 
version. We find that at low energies, 44-48 Gauss quadrature points (24 points for 
k'[ = — Ikyii^n interval and 20-24 points for k" = 2k v ii^i — oo interval and similarly for 
k'l are needed for the discretization of the Kernel of the LS equation, to achieve numerical 
convergence up to 4-decimal places. 

Results and Discussions : We begin the discussions of our results with the singlet 
channel cross sections. The effect of the Ps+H— > e + +H~ channel would be revealed in 
this singlet channel, as H~ belongs to a spin singlet state only. Also, the binding and 
resonance of PsH occur in this channel. In figure-1, we plot the S-wave singlet cross section 
from 0-10 eV to exhibit the influence of this channel (channel-2) over the static-exchange 
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prediction (using both exact and model exchange). In both the cases, we find that this 
channel reduces the singlet scattering cross section significantly, at low energies. This 
implies that this channel effectively makes the attractive potential more strong and we 
can expect a smaller scattering length and a greater binding energy. Using an effective- 
range expansion of the form kcotS = — l/a + r Q k 2 /2 + Dk 4 , and finding solution of the 
equation kcot5 — ik = for the bound state, where k is the momentum of Ps, 5 is the 
S-wave singlet phase shift, r is the corresponding effective range, and D is the coefficient 
for the k 4 term, we obtain the binding energy and scattering length which are tabulated 
in table 1. 

Table 1: Scattering length (S.L.) and binding energies (B.E.) in the singlet channel of Ps-H, 
from static-exchange (SE), static-exchange plus charge-transfer rearrangement (2CH) employing 
exact exchange. Model exchange results without any parameterization (C = 1) are represented 
by (ME), a) present predictions; b) prediction of Campbell et al (H| 





SE 


2CH 


SE(ME) 


2CH(ME) 


22-state 


S.L. 


a) 7.273 

b) 7.25 


6.90 


8.427 


5.544 


5.20 


B.E. 


a) -0.253 

b) -0.263 


-0.291 


-0.139 


-0.394 


-0.634 



As expected, we find the scattering length (S.L.) is reduced and the binding energy 
(B.E.) is increased with the influence of the H~ formation channel. The effect is much 
pronounced in the case of model exchange calculation, thus signifying a lesser role for 
the parameter in the model potential. The scattering length prediction in the 2CH(ME) 
is quite close to the 22-state prediction. Keeping in mind the simple form of H~ wave 
function used in this calculation, and the proven role of the direct r 12 term in the vari- 
ational prediction ||28|| , we get the impression that the impact of this channel in the ab 



initio model could be improved significantly by using a more shopisticated wave function 
for H~ which contains the ri2 term. The small difference in the scattering length (S.L.) 
and binding energy (B.E.) of SE calculations of set (a) and set (b) is due to the fact that 
in (b) we have used an effective range expansion to extrapolate the phase shifts to the 
negative energy region while in set (a) those parameters are directly solved. The positive 
energy scattering phase shifts for (a) and (b) agree exactly. We tabulate the phase shifts 
for different models in table 2 for future reference. 

Table 2: Variation of singlet and triplet scattering phase shifts for different models with 
energy (E = 6.8/c 2 eV): SE-static exchange; 2CH-static exchange plus e + H~ channel; (ME)- 
represents the calculations with model exchange but without any parameterization (C = 1). 



singlet 


triplet 


k 


SE 


2CH 


SE(ME) 


2CH(ME) 


SE 


SE(ME) 


0.1 


2.455 


2.485 


2.387 


2.610 


-0.247 


-0.145 


0.2 


1.927 


1.966 


1.910 


2.176 


-0.489 


-0.283 


0.3 


1.539 


1.577 


1.604 


1.848 


-0.721 


-0.410 


0.4 


1.239 


1.275 


1.381 


1.596 


-0.940 


-0.521 


0.5 


0.997 


1.032 


1.207 


1.396 


-1.143 


-0.613 


0.6 


0.797 


0.831 


1.065 


1.232 


-1.330 


-0.683 


0.7 


0.631 


0.664 


0.946 


1.094 


-1.499 


-0.731 


0.8 


0.491 


0.523 


0.845 


0.977 


-1.653 


-0.757 


0.9 


0.373 


0.406 


0.758 


0.876 


-1.790 


-0.763 


1.0 


0.274 


0.307 


0.683 


0.786 


-1.913 


-0.781 



First we discuss about SE phase shifts where recently a dispute has been raised Table 
2 shows that the present phase shifts for the SE calculation for the singlet and triplet scattering 
agree exactly with the existing SE predictions of Hara and Fraser and Sinha et al fl3| . 
Campbell et al |l(| provided SE phase shifts for other A;— values (Ps energy, E = 6.8k 2 ). At 
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k 2 = 0.1639,0.2478, and 0.5588 (a.u) the singlet phase shifts of Campbell et al are given by 
(table 1 of ref. @) 1.23, 1.00, 0.56, respectively and our predicted values are 1.228, 1.002, and 
0.5614, respectively (all cross sections are in units of 7rag). So all these SE predictions, which use 
BO exchange amplitude as the input, agree among themselves. This suggests that the minimum 
in the BO cross section is a true feature of the model and the results of Adhikari and Mondal 
(AM) []38|| cannot agree with Campbell et al [flOfl . Since, at k 2 = 0.5588, the BO cross sections 
of AM (3^] is about 23.5-7raQ, which is approximately 22% higher than the BO cross sections 
(19.37raQ) of Sinha et al and present calculation. So, it is certain that any static exchange 
prediction made using the BO exchange amplitude of AM as input to the CC equations, will 
disagree with all existing SE results including those of Campbell et al [10]. While trying to get 
an answer to the minimum in the elastic BO cross section, we find that for forward scattering 
the exchange potential changes its sign at about 26-27 eV, which is causing the minimum in the 
cross section. We now come to the main business of the asssement of the influence of the H~ 
formation channel over the SE predictions. 

In table 2, we compare the phase shifts of SE and 2CH calculations. For 2CH we do not 
tabulate the triplet phase shifts since they do not change by the influence of channel-2 (as H~ 
exists only in singlet state). From table 2, we see that singlet phase shifts increase significantly 
with the inclusion of the H~ formation channel. This signifies an increase in the attractive 
potential, which was expected and needed in the CC theory to improve the convergence and 
hence improve the PsH binding energy. It is interesting to note that the effect of this channel 
continues significantly around k = 1.0. 

In figure 2, we plot the elastic cross sections for SE, SE(ME), 2CH, and 2CH(ME) models. 
In this figure we also plot the recent variational prediction |39[ | on the zero energy cross section. 
Elastic cross sections for both the ab initio and model calculations are reduced by the influence 
of the charge transfer channel. At low energies, the 2CH(ME) model gives much lower cross 



section than the 2CH model and the former is quite close to the recent variational prediction [39]. 
We have earlier find that this model can lead to more converged results if the model exchange 
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potential is tuned by means of a parameter C [17, 21]. While trying to analyze the physics 
behind such agreement with measurements and accurate variational predictions, we find that 
effectively, while making approximate mapping of a basis set (VVt( r 2)) belonging to a different 
Fock-space to its original Fock-space (Vv( r i))> basically we are indulging to approximations like: 

^ M (r 2 ) 1 dv 2 = -^ M (ri) (12) 

l r l — r 2| 

(x-r,^-!-^ = ^(x-rx) (13) 



where k^ and k^ are the average values of the square of the momenta of the respective electrons 
when they were bound in the hydrogen (?/^(r 2 )) or in th.6 posit roniiim atom. {(j)y (x - ri)) and 
C is a parameter which is fixed to unity but could be varied to tune results. Clearly, the 
above relations are introducing some amount of e\ — e 2 correlation and the continuum effects 
of the target (eqn.()) and the projectile (eqn.()), respectively. Also, these effects would be 
modified when the value of CkZ or Ckf, are changed. Interestingly, the model provides best 



results [17, 22, 15] when Ck^ and Ck„ become close to the ionization energies for the respective 
electrons. Although the mixing of such model exchange with ab initio charge-transfer process 
yielded very good response, we are doubtful about the overcompleteness of the Hilbert space if 
such mixing is continued with a larger basis for Ps. However, from the results of the 2CH model 
and from our experience in such rearrangement for positron-atom scattering [], we understand 
that for the full ab initio model (2CH) the possibility of overcompleteness is remote unless 
a very large basis is used. The information we gain here in terms of the convergence of the 
2CH prediction appears consistent and encouraging due to the fact that virtual excitations of 
Ps are absent here. The charge-transfer reaction becomes exothermic for n > 5 and higher 
discrete excited states and continuum of Ps and hence the effect of such a channel is likely to 
be significant with higher Ps eigen states and pseudostates. Consideration of this channel along 
with the Ps-excited states would also reveal its effect on resonances, which manifest in the theory 
with the inclusion of Ps excited states. 

Another interesting feature of such study is that it provides the H~ formation cross section 
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which is of importance in Astrophysics. Here we provide the H formation cross sections from 



Ps(ls)-H(ls) scattering in figure-2. Solid curve is obtained with the regularized exchange [21, 18| 
while the dotted curve is obtained with ab initio exchange. It is observed that from 60 eV 
onwards the two sets of results are same signifying the fact that the regularized model exchange 
asymptotically coalesce with the exact one. From the two curves we also see that, in the low 
and intermediate energies, the cross sections calculated with the regularized exchange CC model 
(2CH(ME)) are higher than those calculated with ab initio exchange (2CH). This was quite 
expected as the regularized exchange potential, as discussued above, was found to simulate the 
continuum effect substantially and yield converged results for the PsH binding and resonances 
1 17, 22 1 . The contribution from virtual Ps excitations, no doubt will increase the formation 
cross section and this trend has been revealed in the 2CH(ME) calculation. 

Conclusion : In summary, we make a first time study of the effect of charge transfer recom- 
bination for Ps scattering from hydrogen. Considering the rearrangement of PsH to e + H~ we 
get the information that such rearrangement can introduce the continuum effect more effectively 
than considering only exchange. The attractive potential introduced by this channel improves 
the scattering length and binding energy predictions significantly. The influence of the channel 
is also found to help converging the low energy cross sections (see figure 1). Similar effect of 
charge-transfer recombination was found in the cases of positron-atom scattering [jl], || and in 
the present case of Ps-atom scattering it could be more interesting as the threshold for the 
transfer process Ps(ls)+H(ls)— > e + +H~ decreases with higher Ps states and the reaction be- 
comes exothermic for n > 5 discrete excited states and continuum of Ps. So, it is expected that 
the channel would contribute substantially through virtual Ps excitations and lead to further 
convergence for the low energy Ps scattering parameters. As a by product we get H~ formation 
cross sections which has astrophysical importance. 

The work is supported by the Fundacao de Amparo a Pesquisa do Estado de Sao Paulo of 
Brazil via project number 99/09294-8. 
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Appendix I 



The Born-Oppenheimer exchange matrix element for the elastic transition is given by: 



iMrL+x)^^ _ x)(j)ls ( r2 ) dxdridr2 . (14) 



lsls.lsls 

^/i J J J 

x e"2 



where kj, are the initial and final momentum of the Ps atom with respect to the center of 
mass. We use Ti = c = mi = 1, m 2 is the mass of the electron or positron; the Ps mass m\ = 2. 
Splitting with the four terms of V^t of eqn(4), we write Bf sls lsls = B\ + Bi + -B3 + -B4, where 
subscripts denote the terms on the r.h.s. of eqn(4), in sequence. We use Fourier transform, 
delta-function integration, and Chasire integrals to reduce each nine dimensional integrals to 
two dimensions. We follow same procedure for Bi,Bi,B^, B4. We detail for B4, which contains 
l/|ri — rz\ term of V^J. Using Xi s (r) = -j=e~P r and 0i s (r) = -^=e~ ar , with /3 = 0.5 and 



where 



a = 1.0, we perform integration over dx first, then over dr±, and finally over dr? to get, 

5 4 = 4/? 2 a fduuil-u) [\ W (±°) 2 1-(1-°)1- / 2 « 2 (15) 
Vo Jo Kill dps. J III V/U 2 <9/i 2 7 (|A| 2 +Hla) 

£ = /3 2 + u (l-u)|0.5(ki-k f )| 2 (16) 

/1 2 . = va 2 + (l--u)// 2 + -u(l--u)|Q| 2 (17) 

Q = (1 - 0.5u)ki - 0.5(1 - u)kf (18) 

A = (0.5u + w - 0.5«t;)ki - 0.5(1 + u + v - uv)\a t (19) 

M2a = M2 + a (20) 

For the rest of -Bi,-B 2 and B3, we perform integration over dr\ first, then over dx, and finally 
over dr2 and arrive at similar final forms as above, with changed definition for ^2, Q, A etc. 
Integration over each of du and dv is converged with mere eight Gauss-quadrature points. 
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Appendix II 



The charge transfer matrix element for the Ps+H^ e + +H channel is given by: 



Bf s2;lsls (k { M) = -^ J J /e- k /- x V^-(ri,r 2 )^ 

where kj, now represents the momentum of the ejected positron with respect to the center 
of mass rest on the proton; mi = 1; ipH-( T ii r 2) = jj(e~ aTl e~ br2 + e~ aT2 e~ bri ^j , the H~ wave 
function with binding energy E = —0.51330 and a = 1.03925, b = 0.28309 and the normalization 
constant has been worked out to be N = 31.80348105. As in Appendix I, we split the integral 
as Bf g2 lsls = B\ + B2 + S3, with the three terms of V^} of eqn(5), where subscripts denote 
the terms on the r.h.s. of eqn(5), in sequence. Using the symmetry of the H~ wave function we 
write Bj = Bj(a,b) + Bj(b, a); j = 1,2,3. We use Fourier transform, delta-function integration 
and Chasire integrals and evaluate B2 exactly. 

B , >\ 64y/27r 2 a 

2l °' ; " (« + 6) 3 (|qi/ 2 | 2 + /? 2 ) 2 + (|q| 2 + a 2 ) 2 1 ' 

where qi/2 = 5^ — kf and q = ki — kf . Integration over dr2 in B\ is simple and in S3 (a, b), it 

gives 



/ 



e -(a+6)r 2 dr 2 = 87T _ 87re-^+ b > _ 4^+^ 

|x-r 2 | (a + b) 3 x (a + b) 3 x (a + b) 2 1 ' 



With this form of /, we write S 3 (a, b) = S3 (a, b) + Sf(a, b) + Sf (a, 6). Now S3 (a, 6) cancels out 
with Bi(a,b) and similarly for B\(b,a) and Si (6, a). We evaluate S 2 (a, 6) and obtain Sf(a, 6) 
from it with a derivative w.r.to (a + 6). The final form of S 2 (a, 6) is 

B»(a,b) = — ■; p-V / dun -— - fl ^ lo . r^pr (24) 

(a + 6) 3 Jo V/^<W /i{|Q| 2 + (/U + a) 2 } 2 v ; 

M 2 = u p 2 + (l- u )( a + b) 2 +u(l-u)\q_ 1/2 \ 2 (25) 
Q = (1 - 0.5n)ki - (1 - u)k f (26) 



where 
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Figure Captions: 

Fig. 1: Variation of S-wave singlet cross sections (in units of vtoq) for Ps-H scattering 
employing static-exchange (SE), SE plus charge-transfer rearrangement to e + H _ (2CH). (ME) 
represents calculations using model exchange plf : 

Fig. 2: Variation of elastic cross sections (in units of ttoq) for Ps-H scattering employing 
static-exchange (SE), SE plus charge-transfer rearrangement to e + H~ (2CH). (ME) represents 
calculations using model exchange ||21|| : 

Fig. 3: H formation cross sections from Ps(ls)-H(ls) scattering (in units of vtoq)- 2CH 
corresponds to results fron the two-channel calculation with ab initio exchange while 2CH(ME) 
is the same using regularized model exchange. 
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